OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



Characteristics of HIV-1 Natural Drug Resistance- 
Associated Mutations in Former Paid Blood Donors in 
Henan Province, China 

Yi Bao, Di Tian, Ying-Ying Zheng, Hong-Li Xi, Dan Liu, Min Yu, Xiao-Yuan Xu* 

Department of Infectious Diseases, Peking University First Hospital, Beijing, China 



Abstract 

Background: Natural drug resistance is a major cause of antiviral treatment failure. The characteristics of HIV-1 natural drug 
resistance-associated mutations in former paid blood donors in Henan Province remain unclear. 

Methods:One hundred and fifty HIV-1 -positive plasma samples were collected. Plasma viral RNA was extracted for pol gene 
amplification and sequencing. The sequencing results were submitted to the HIV-1 drug resistance database for drug- 
resistance analysis. 

Results: The rates of natural drug resistance and resistance-associated mutations were 17.7% (19/107) and 40.2% (43/107), 
respectively. The rates of PI major, PI minor, NRTI, and NNRTI mutations were: 0, 30.8% (33/107), 10.3% (1 1/107), and 18.7% 
(20/107), respectively. Nine cases (8.4%) had both NRTI and NNRTI resistance-associated mutations. Seven cases (6.5%) had 
PI minor, NRTI and NNRTI resistance-associated mutations. NNRTI resistance was the most serious, followed by NRTI 
resistance and PI resistance. Polymorphism mutation sites with mutation rates in the protease region higher than 60.0% 
were: L63A/P/S/T 89.7%, V77I 82.2%, I72E/M/K/T/V 80.4%, I93L 75.7%, and E35D 72.9%. Polymorphism mutation sites with 
mutation rates in the RT region higher than 60.0% were: I135A/L/M/R/T/V 93.5%, T200A/E/I/P/V 89.7%, Q278E/K/N/T 88.8%, 
S162C/Y 82.2%, and K277R/S 66.4%. The distribution of 107 gene sequences was scattered, with some drug-resistant strains 
grouped in the same cluster. 

Conclusion: The natural drug resistance mutation rate of HIV-1 in former paid blood donors in Henan Province was 17.7%, 
with NNRTI resistance the most serious. The distribution of drug-resistant strains was scattered, with some correlations 
found in certain resistance loci. 
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Introduction 

Under natural circumstances, HIV can destroy the human 
immune system 10—15 years after the development of AIDS. At 
this stage, the patients may display a variety of opportunistic 
infections, tumors or other serious clinical symptoms [1—3], Up to 
now, there are no effective drugs and therapies for complete cure 
of AIDS. Extensive application of ART has greatly reduced the 
morbidity and mortality of HIV/ AIDS patients [4-6]. 

According to the AIDS epidemic report from Joint United 
Nations Programme on HIV/AIDS(UNAIDS) and the World 
Health Organization (WHO) [7], there were 34 million HIV 
carriers and AIDS patients, 2.7 million new HIV infections and 
1 .8 million deaths worldwide in 20 1 0. By the end of 20 1 1 , in low- 
and middle-income countries, a total of more than 8 million HIV 
carriers received antiretroviral therapy, of which about 562,000 
were children. During 2003-2011, the number of patients 
receiving antiretroviral therapy increased by 20 times in develop- 



ing countries. During 2010-201 1, the number of patients needing 
antiretroviral treatment increased by 20%. Based on the 
mechanism of action, antiviral drugs are divided into different 
groups, which include protease inhibitors, reverse transcriptase 
inhibitors, integrase inhibitors, fusion inhibitors, entry inhibitors 
and mixture of different agents. Among them, protease inhibitors 
and reverse transcriptase inhibitors are used more often in 
clinics [8-11]. 

However, HIV genome replication is error-prone. The virus 
exists in patients as a large and complex pool of quasispecies form 
[12,13]. The selection pressure of antiviral drugs causes the 
accumulation of resistant strains that gradually become the 
dominant strains, leading to resistance to protease inhibitors and 
reverse transcriptase inhibitors. In addition, a resistance mutation 
caused by one drug may confer resistance to another drug that is 
not used, a phenomenon called cross-resistance [14- 17]. Numer- 
ous studies have shown that drug resistance mutation is one of the 
main reasons leading to clinical antiviral treatment failure. 
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Furthermore studies in the United States and Western Europe 
showed that the spread of drug-resistant strains increased 
significandy, and about 10.7% ~ 27.6% of new infections were 
caused by drug-resistant strains of HIV-1 [18— 19]. Natural drug 
resistance is a major cause of antiviral treatment failure. Before 
antiviral treatment, drug resistance tests for HIV-1 strains can 
guide the clinical medication, which will greatly reduce the risk of 
initial treatment failure [20-25]. 

Since 2003, most patients undergoing AIDS antiretroviral 
treatment were former paid blood donors in Henan Province. 
HIV drug resistance caused by antiretroviral therapy has become 
a major obstacle to clinical treatment. In recent years, drug 
resistance data in Henan Province have been reported [2 6-2 7], but 
the research about primary drug resistance is lacking. With 
infection rates rapidly rising, the primary drug resistance data also 
need to be constantly updated and supplemented. 

In this study, we investigated the characteristics of HIV-1 
natural resistance mutations in former paid blood donors in 
Henan Province. Understanding primary drug resistance before 
antiretroviral treatment may facilitate the development of effective 
treatment programs in order to achieve the greatest degree of viral 
suppression, and understand molecular evolution traceability rules 
of HIV-1 resistance genes. It is also important for reducing the 
prevalence of drug-resistant strains and the effective use of drugs. 

Methods 

Ethics Statement 

This study was approved by the Ethical Committee of Beijing 
You'an Hospital Affiliated to Capital Medical University and 
written informed consent was obtained from all participants. Each 
consented participant was given a unique number. The biological 
and behavioural information was linked anonymously to protect 
the participants' privacy. This procedure was approved by the 
ethics committee. 

Samples 

Between 2005 and 2007, 150 plasma samples were collected in 
Beijing YouAn Hospital from paid blood donors in Henan 
Province. All patients were confirmed as HIV-1 positive using 
Western blot. The ages of patients were between 28 and 60 years, 
with a mean age of 44.7 years. There were 86 male cases, 
accounting for 57.3%, and 64 female cases, accounting for 42.7%. 
HIV infected peripheral blood samples were collected using 
EDTA anticoagulant tubes before treatment with antiviral drugs. 
Plasma samples were stored at — 80°C. 

Nucleic Acid Extraction 

HIV-1 RNA was extracted from 200 ul of plasma using the 
MiniBEST Viral RNA/DNA Extraction Kit (TaKaRa, Dalian, 
China) according to the manufacturer's instructions. 

Nested Polymerase Chain Reaction 

HIV-1 Pol gene was amplified using nested-PCR. The TaKaRa 
PrimeScript One Step RNA RT-PCR Kit(TaKaRa, Dalian, 
China) was used for the first round of PCR. TaKaRa Ex Taq 
Kit(TaKaRa, Dalian, China) was used for the second round of 
PCR. The primers used and the experimental procedure were 
described previously [28,29]. The primer sequences are shown in 
Table 1. The tests used disposable supplies (pipette tips, gloves, 
etc), and were carried out by skilled personnel. Two negative 
controls were included for each experiment. The positive results 
were repeated twice. 



Electrophoresis and Sequencing 

For electrophoresis, 5 ul of first round PCR product was loaded 
onto a 2.0% agarose gel. The gel was run with a constant voltage 
of 120V for 25min. PCR positive products were used for 
sequencing analysis(SinoGenoMax Co, Ltd, Beijing). The ampli- 
fication products were 1315-bp in length. All the gene sequences 
were uploaded to the Genebank. 

HIV-1 Natural Drug Resistance-Associated Mutation 
Analysis 

MEGA4.02 software was used to analyze the drug resistance 
gene sequences. The phylogenetic tree was constructed using the 
Bootstrap Neighbor Joining method. Genotyping was conducted 
using BLAST tool provided by the U.S. Los Alamos HIV database 
(http://www.hiv.lanl.gov). The edited sequences were submitted 
to the Stanford HIV Drug Resistance Database (http://hivdb. 
stanford.edu/index.html) (Version 6.2.0, last updated 05/29/12) 
for analysis of drug resistance, resistance loci and resistance- 
associated mutations. The preliminary target sequences were 
evaluated for drug sensitivity, potentially low resistance, low 
resistance, medium resistance or high resistance. 

Results 

General 

One hundred and seven HIV-1 pol gene sequences were 
obtained. The positive rate was 71.3% (107/150). The ages of 
patients were between 29 and 60, with an average age of 45.0. 
Among them, 66 (61.7%) were males, with an average age of 45.4 
years, and 41 (38.3%) were females, with an average age of 44.6 
years. As determined by Stanford HIV Drug Resistance Database, 
40.2% (43/107) samples had one or more resistance-associated 
mutations. These patients had an average age of 47.6 years, with 
male accounting for 69.8% (30/43, mean age 47.7 years), and 
female accounting for 30.2% (13/43, mean age 47.5 years). The 
CD4 cell counts, plasma HIV RNA levels and HCV RNA levels 
were determined, and the results were summarized in the table 2. 
However, there was not enough blood plasma for HBV 
measurement. 

Phylogenetic Analysis of the Characteristics of HIV-1 
Natural Resistance Mutations 

As shown in Figure 1, the distribution of 107 sequences was 
scattered. They were mainly clustered together with reference 
strains B.CN.2001.CNHN24.AY180905 and B.FR 1983.HXB2- 
LAI-IIIB-BRU.K03455. The subtype of these strains was B 
subtype. Samples with black spots represented the drug-resistant 
mutations. HNHIV3, HNHIV50, and HNHIV139 were clustered 
together, and had LI 01 mutation. HNHIV46, HNHIV74, 
HNHIV77, and HNHIV137 were clustered together, with 
HNHIV46 and HNHIV137 having A71V, D67N, K101E, and 
G190A mutations, HNHIV74 having A71T mutation, and 
HNHIV77 having A71V and K103N mutations. Samples 
HNHIV6, HNHIV58, and HNHIV149 were clustered together, 
and had A71T mutation. Samples HNHIV19, HNHIV68, 
HNHIV69, HNHIV93, HNHIV97, HNHIV113, and 
HNHIV114 were clustered together, with HNHIV19 and 
HNHIV93 grouped in a smaller cluster. The drug resistance 
mutation loci included A71V, D67H, T69N, K70R, T215F, 
K219E, A98G, V179D, and Y181C, with HNHIV97 and 
HNHIV113 clustered together. HNHIV97 locus mutations 
included D67N, T69A, K70R, T215Y, K219E, K103N, and 
Y181C. HNHIV113 site mutations included A71T, M41L, 



PLOS ONE | www.plosone.org 



2 



February 2014 | Volume 9 | Issue 2 | e89291 



Characteristics of HIV-1 Drug Resistance 



Table 1. Primers used in this study. 





Nested PCR 


Name 


Sequences(5 -3 ) 


Position (HXB2) 


1st PCR 


MAW-26 


TGGAAATGTGGAAAGGAAGGAC 


2147-2166 




RT-21 


CTGTATTTCTGCTATTAAGTCTTTTGATGGG 


2539-2519 R 


2nd PCR 


PRO-1 


CAGAGCCAACAGCCCCACCA 


2946-2961 




RT-20 


CTGCCAGTTCTAGCTCTGCTTC 


2050-2029 R 


For sequencing 


PRO-1 


CAGAGCCAACAGCCCCACCA 


3539-3509 




RTA 


GTTGACTCAGATTGGTTGCAC 


2147-2166 




RTB 


CCTAGTATAAACAATGAGACAC 


3462-3441 



doi:1 0.1 371 /joumal.pone.0089291 .t001 



M184V, L210W, T215Y, K103N, and Y181C. Samples 
HNHIV68, HNHIV69, and HNHIV1 14 were clustered together, 
with both HNHIV68 and HNHIV69 having mutations K103N 
and Y181C, HNHIV114 having mutation A71V. Samples 
HNHIV26, HNHIV56, and HNHIV145 were clustered together, 
with HNHIV26 having mutations L10I, A71T, K103N, and 
Y181C, and HNHIV56 and HNHIV145 having mutation A71T. 
Samples HNHIV4, HNHXV83, HNHIV100, and HNHIV121 
were clustered together, and had A71T mutation. In addition, the 
PCR sequences with < 1 % difference are 5 1 and 1 40, 56 and 145. 

In samples with LI 01 mutation, HNHIV24 was not clustered 
together with HNHIV3, HNHIV50 and HNHIV139. In samples 
with A71T mutation, HNHIV6, HNHIV58, and HNHIV149 
were clustered together, HNHIV4, HNHIV83, HNHIV100, and 
HNHIV121 were clustered together, and HNHIV10 and 
HNHIV74 had scattered distribution. In samples with A71V 



mutation, HNHIV23, HNHIV44, HNHIV63, and HNHIV114 
had scattered distribution. All samples with both A71V and 
K103N mutations had scattered distribution. In samples with the 
same mutation, HNHIV19 and HNHIV93 were clustered 
together, and HNHIV65 was not clustered together with 
HNHIV68 and HNHIV69. 

The Distribution of Natural Resistance Locus 

PI had no major resistance-associated mutations. PI minor 
resistance-associated mutation rate was 30.8% (33/107). The 
mutation positions were A71T, L10I, A71V, and Q58E, with 
A71T and LI 01 happening simultaneously, leading to NFV 
potential low resistance. Q58E led to TPV/r potential low-level 
resistance, with a resistance rate of 1.9% (2/107). 

NRTI resistance-associated mutations were found in 1 0 patients 
(9.3%). Nine were resistant to NRTI, with a resistance rate of 



Table 2. Characteristics of treatment-naive HIV-infected study population in former paid blood donors in Henan province, China. 



Characteristics 


Total n (%) 


HIV-1 pol gene sequences 
were Obtained n (%) 


HIV drug resistance n 

(%) 


Sex 


150 


107 


43 


Male 


86(57.3) 


66(61 .7) 


30(69.8) 


Female 


64(42.7) 


41 (38.3) 


13(30.2) 


Median age at inclusion, years (range) 


44.7 (29-60) 


45.0(29-60) 


47.6(30-60) 


Baseline CD4 + T cell count (cells/ml) 


73 


54 


21 


>500 


7(9.6) 


5(9.3) 


0 


350-499 


38(52.1) 


26(48.1) 


12(57.1) 


200-349 


28(38.3) 


23(42.6) 


9(42.9) 


50-199 


0 


0 


0 


<50 


0 


0 


0 


HIV Baseline viral load (copies/ml plasma) 


150 


107 


43 


<1000 


2(1.3) 


0 


0 


1 000-9999 


26(17.3) 


3(2.8) 


0 


1 0 000-99 999 


85(56.7) 


67(62.6) 


27(62.8) 


>100 000 


37(24.7) 


37(34.6) 


16(37.2) 


HCV Baseline viral load (copies/ml plasma) 


150 


107 


43 


<1000 


54(36.0) 


39(36.5) 


13(30.2) 


1 000-9999 


3(2.0) 


0 


0 


1 0 000-99 999 


9(6.0) 


3(2.8) 


2(4.7) 


>100 000 


84(56.0) 


65(60.7) 


28(65.1) 



doi:1 0.1 371 /joumal.pone.0089291 .t002 
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Figure 1. Phylogenetic analysis the characteristic of HIV-1 
natural resistance mutations in former paid blood donors in 
Henan Province. The bootstrap values were performed by 5000, 
greater than 80 clusters shown in figure, black spots represent the drug- 
resistant mutations samples, black triangles represent the aggregation 
of the same cluster, reference strains for subtype B are: 
B.CN.2001 .CNHN24.AY1 80905 and B.FR.1 983.HXB2-LAI-IIIB-BRU.K03455. 
doi:10.1371/journal.pone.0089291.g001 

8.4% (9/107). In these cases, 4.7% (5/107) were high resistance, 

I. 9% (2/107) were intermediate resistance, 6.5% (7/107) were low 
resistance and 1.9% (2/107) were potential low resistance. Single 
point mutation Ml 84V led to high resistance to 3TC and FTC, 
low drug resistance to ABC, and potential low resistance to DDI; 
T2 1 5Y caused intermediate resistance to AZT and D4T, and low 
resistance to ABC, DDI and TDF; M184L did not cause resistance 
to NRTI. When comparing HNHIV97 and HNHIV93, T215F 
was replaced by T2 1 5Y, which changed intermediate resistance to 
DDI to low resistance. 

NNRTI resistance-associated mutations were found in 20 cases 
(18.7%). Eighteen cases were NRTI resistant, with a resistant rate 
of 16.8% (18/107). Among them, 15.9% (17/107) were high 
resistance, 10.3% (11/107) were intermediate resistance, 1.9% (2/ 
107) were low resistance, and 0.9% (1/107) were potential low 
resistance. In seven cases Y181C occurred together with K103N, 
causing high resistance to EFV and NVP, and intermediate 
resistance to ETR and RPV. The mutations of both K103N and 
VI 081 caused high resistance to EFV and NVP. Y181C caused 
intermediate resistance to ETR and RPV. K101E and G190A 
caused high resistance to NVP, intermediate resistance to EFV, 
and low resistance to ETR and RPV. V179E did not cause 
resistance to NNRTI. Seven patients (6.5%) had PI minor, NRTI 
and NNRTI resistant mutants. Nine patients (8.4%) had both 
NKTI and NNRTI resistance mutations. Specimens HNHIV19 
and HNHIV93 had mutations at the same locus, and HNHIV46 
and HNHIV137 had mutations at the same locus. The results 
were summarized in the table 3. 

HIV Drug Resistance Mutation Annotations 

In PI minor resistance-associated mutations, LI 01 accounted for 
4.7%, which is slightly lower than the rate (5%— 10%) of untreated 
population in the Stanford HIV drug resistance database 
annotation. A71T and A71V mutation rates were 14.9% and 

II. 2%, respectively. These were polymorphic mutations, and 
occurred at a higher rate than the rate (2—3%) of the untreated 
population in the database annotation. 

In the RT region, among mutations associated with NRTI 
resistance, M41L often occurred together with T215Y, which 
caused high resistance to AZT and d4T, and low resistance to 
DDI, ABC and TDF. D67N accounted for 2.8% and caused AZT 
and d4T resistance. T69A and T69N accounted for 0.9% and 
1.9%, respectively. Mutations at this site are usually caused by 
drug selection pressure. But the samples were collected before 
antiviral treatment in the present study. M184L was a rare 
mutation, which accounted for 0.9%. M184V accounted for 1.9%, 
and was responsible for high resistance to 3TC and FTC, and low- 
level resistance to ddl and ABC. L210W was often associated with 
M41L and T215Y mutations, and caused all NRTIs except 3TC 
and FTC resistance. K219E accounted for 2.8%, and could 
decrease AZT sensitivity in the presence of K70R or T215Y/F. 
Among mutations associated with NNRTI resistance, A98G 
accounted for 1.9%, and decreased the NVP and EFV sensitivity. 
K103N accounted for 12.1%, and caused high-level resistance to 
NVP and EFV. V108I accounted for 1.9%, and reduced NVP and 
EFV susceptibility. V179D/E reduced NVP and EFV suscepti- 
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Table 3. HIV-1 drug resistance mutations natural distribution of former paid blood donation population in Henan province. 



Number 


Gender 


Age 


Drug Resistance Interpretation: 
PR 


Drug Resistance Interpretation: RT 












PI Maj SpeR PI Min 


Spe R 


NRTI 


Spe R 


NNRTI 


Spe R 


HNHIV1 


female 


53 


A71T 












HNHIV3 


male 


51 


L10I 












HNHIV4 


female 


46 


A71T 












HNHIV6 


female 


36 


A71T 








_ 


_ 


HNHIV10 


male 


50 


A71T 


- 


- 


- 


_ 


_ 


HNHIV13 


female 


37 


A71V 




Ml 84V 


H:3TC,FTC LABC P:DDI 


K103N, 
Y181C 


H:EFV,NVP l:ETR,RPV 


HNHIV15 


male 


57 






K70N 


P:3TC,ABC,D4T,DDI,FTC,TDF 


VI 081 


P:NVP 


HNHIV19 


female 


49 


A71V 




D67H, T69N, K70R, 
T215F, K219E 


H:AZT,D4T l:DDI,TDF LABC 


A98G, 

V179D, 

Y181C 


H:NVP l:EFV,ETR,RPV 


HNHIV20 


male 


56 


A71V 


_ 


_ 


_ 


K103N 


H:EFV,NVP 


HNHIV23 


male 


46 


A71V 












HNHIV24 


male 


54 


L10I 


- 


- 


- 


- 


- 


HNHIV26 


male 


49 


L10I, 
A71T 


P:NFV 


- 


- 


K103N, 
Y181C 


H:EFV,NVP l:ETR,RPV 


HNHIV29 


male 


37 


A71T 




T215Y 


I:AZT,D4T L:ABC,DDI,TDF 


K103N 


H:EFV,NVP 


HNHIV36 


male 


51 


Q58E 


P:TPV/r 










HNHIV39 


male 


49 


A71V 


- 


- 


- 


K103N 


H:EFV,NVP 


HNHIV44 


male 


51 


A71V 












HNHIV46 


male 


47 


A71V 




D67N 


LAZT.D4T 


K101E, 
G190A 


H:NVP l:EFV L:ETR,RPV 


HNHIV50 


male 


41 


L10I 








- 


- 


HNHIV56 


male 


49 


A71T 




- 


- 






HNHIV58 


male 


46 


A71T 












HNHIV63 


male 


46 


A71V 








- 




HNHIVdd 


female 


47 










K103N, 
Y181C 


■ ■ r- r-i i iii ip, i.rrn nni; 

H:brv,NVr l±IK,KPV 


HNHIV68 


male 


52 










K103N, 
Y181C 


H:EFV,NVP l:ETR,RPV 


Number 


Gender 


Age 


Drug Resistance Interpretation: PR 


Drug Resistance 
Interpretation: RT 














PI Maj SpeR PI Min 


Spe R 


NRTI 


Spe R 


NNRTI 


Spe R 


HNHIV69 


male 


53 










K103N, 
Y181C 


H:EFV,NVP l:ETR,RPV 


HNHIV70 


male 


60 






- 


- 


K103N, V108IH:EFV,NVP 


HNHIV74 


male 


49 


A71T 












HNHIV76 


male 


58 










K103N 


H:EFV,NVP 


HNHIV77 


male 


50 


A71V 




- 


- 


K103N 


H:EFV,NVP 


HNHIV79 


female 


54 






M184L 








HNHIV83 


female 


53 


A71T 












HNHIVod 


male 


36 










V179E 




HNHIV93 


female 


35 


A71V 




D67H, T69N, K70R, 
T215F, K219E 


H:AZT,D4T l:DDI,TDF LABC 


A98G, 

V179D, 

Y181C 


H:NVP l:EFV,ETR,RPV 


HNHIV97 


male 


56 






D67N, T69A, K70R, 
T215Y, K219E 


H:AZT,D4T l:TDF LABC.DDI 


K103N, 
Y181C 


H:EFV,NVP l:ETR,RPV 


HNHIV100 


female 


47 


A71T 












HNHIV113 


male 


30 


A71T 




M41L, M184V, 
L210W, T215Y 


H:3TC,ABC,AZT,D4T,DDI,FTC 
l:TDF 


K103N, 
Y181C 


H:EFV,NVP l:ETR,RPV 


HNHIV114 


female 


39 


A71V 
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Table 3. Cont. 







Mi imkar 
IM UIII Ucl 


Gender 


Age 


Drug Resistance Interpretation: 

PR Drug Resistance Interpretation: RT 








PI Maj 


SpeR PI Min Spe R NRTI Spe R NNRTI Spe R 


l imi iii/i i r 

HNHIV1 1 6 


female 


54 




A71T - - - 


HNHIV121 


male 


36 




A71T - 


HNHIV137 


male 


49 




A71V - D67N L:AZT,D4T K101E, H:NVP l:EFV LETR.RPV 

G190A 


HNHIV139 


male 


51 




L10I - - - 


HNHIV145 


male 


40 




A71T - - - - 


HNHIV149 


female 


52 




A71T - - - 


HNHIV150 


male 


48 




- - - V179E 


PI mar, PI Major Resistance Mutations; PI min, PI Minor Resistance Mutations; Spe R, Specific resistance; H, High-level resistance; I, Intermediate resistance; L, Low-level 
resistance; P, Potential low-level resistance; ATV/r, atazanavir/r; DRV/r, darunavir/r; FPV/r, fosamprenavir/r; IDV/r, indinavir/r; LPV/r, lopinavir/r; NFV, nelfinavir; SQV/r, 



saquinavir/r; TPV/r, tipranavir/r; 3TC, lamivudine; ABC, abacavir; AZT, zidovudine; D4T, stavudine; DDI, didanosine; FTC, emtricitabine; TDF, tenofovir; EFV, efavirenz; ETR, 
etravirine; NVP, nevirapine; RPV, rilpivirine; Mutation sites according to the area of the amino acid position in the form of so-called "Ml 84V" records, said sharing 
sequence of RT's 184th "M" was replaced by "V". 
doi:1 0.1 371 /journal.pone.0089291 .t003 



bUity by about 2-fold. K103R and V179D mutations reduced 
NVP and EFV susceptibilities by about 10-fold. Y181C accounted 
for 8.4%, and caused high-level resistance to NVP. G190A 
accounted for 1.9%, and caused high level resistance to NVP and 
intermediate resistance to EFV. The results were summarized in 
the table 4. 

Natural Resistance Rate of HIV-1 in Former Paid Blood 
Donation Population 

In the protease region, NFV and TPV/r potential resistance 
and low resistance rates were both 0.9%. In the RT region, NRTI 
resistance-associated mutations were mainly high resistance, with 
AZT and D4T high resistance rate of 3.7%. The next most often 
mutation was low resistance, with ABC low resistance rate of 
4.7%. High resistance did not occur with TDF, but intermediate 
resistance rate was as high as 3.7%. Intermediate resistance and 
low resistance did not occur with 3TG and FTC. NRTI resistance 
rates were: 3TC 2.8% (3/107), ABC 6.5% (7/107), AZT6.5% (7/ 
107), D4T 7.5% (8/107), DDI 6.5% (7/107), FTC 2.8% (3/107), 
and TDF5.6% (6/107). Overall NNRTI resistance was mainly 
high resistance, with EFV and NVP resistance rates of 10.3% and 
14.0%, respectively. The next most often mutation was interme- 
diate resistance, with ETR and RPV intermediate resistance rates 
both of 8.4%. Low resistance and potential low resistance did not 
occur with EFV. High resistance and potentially low resistance did 
not occur with ETR and RPV. Intermediate resistance and low 
resistance did not occur with NVP. NNRTI resistance rates were: 
EFV 14.0% (15/107), ETR 10.3% (11/107), NVP 15.4% (16/ 
107), and RPV 10.3% (1 1/107) (Table 5). NNRTI resistance was 
the most serious, followed by NRTI resistance. PI resistance was 
the least serious. 

The Characteristics of Other Mutations 

Polymorphic mutations were mainly in the protease region. The 
sites with mutation rates higher than 60.0% included: L63A/P/S/ 
T 89.7% (96/107), among which L63P was 70.1% (75/107); V77I 
82.2% (88/107); I72E/M/K/T/V 80.4% (86/107), among which 
I72V was 63.6% (75/107); I93L 75.7%; E35D 72.9%. In addition, 
there were some other mutation sites: I62V 24.3%, T12A/K/P/S 
19.6%, R41K 15.9%, N37D/S 11.2%, I64L/M/V 10.3%, R57K 
10.3%. 



The reverse transcriptase region also contained polymorphic 
mutations. The sites with mutation rates higher than 60.0% 
included: I135A/L/M/R/T/V 93.5% (100/107), among which 
I135V was 80.4% (86/107); T200A/E/I/P/V 89.7% (96/107), 
among which T200A was 81.3% (87/107); Q278E/K/N/T 
88.8% (95/107), among which Q278E was 85.0% (91/107); 
S162C/Y 82.2% (88/107), among which S162C was 81.3% (87/ 
107); K277R/S 66.4% (7 1/ 107), among which K277R was 64.5% 
(69/107). Other mutation sites included: I293V 45.8%, V245E/ 
G/I/K/M/L/Q/T 28.0%, R211E/G/M/K 26.2%, F214I 
14.0%, K122E 14.0%, A272G/P/S 13.1%, Q207E/H/K 1 1.2%. 

Discussion 

HIV- 1 reverse transcriptase lacks the proof-reading function. As 
a result, it will produce a base mismatch roughly every 10, 000-30, 
000 nucleic acids during the synthesis. Results from the HIV-1 
drug resistance mutation research by the International AIDS 
Society-USA (updated in March 2013) have revealed that PI 
resistance mutation sites are LI 01, K20M, V32I, M36I, M46I/L, 
I47V/A, I50V, Q58E, A71V, G73S, V82A/F/T, I84V, 
L89V,L90M; NRTIs resistance mutations are M41L, A62V, 
K65R, D67N, K70E/R, Y115F, Ml 84V A, L210W, T215Y/F, 
K219Q/E; and nNRTIs resistance mutations are V90I, LI 001, 
K101E/P, K103N/S, V106A/I/M, E138A/G, Y181C/I/V, 
Y188C/L, G190S/A, M230L [24]. 

In Henan Province, one of the earliest areas in China with free 
anti-AIDS treatment, drug resistance caused by long-term 
treatment has long been a big concern. Drug-resistant strains of 
HIV in the region have certain characteristics. This study is a 
cross-sectional survey carried out on the drug-resistant strains of 
HIV in Henan Province, there are 43 drug resistance patients is 
unknown, the actual rate of primary drug resistance all 150 
patients may also differ. All of 107 HIV-1 pol gene sequences were 
subtype B. This is consistent with the fact that subtype B is the 
main strain in patients infected with HIV through blood 
transmission in rural areas of central China[30-33]. 

In the phylogenetic tree, the distribution of 107 gene sequences 
was scattered, with some drug-resistance strains grouped in the 
same cluster. But in the same cluster, resistance loci were not 
identical, with the same resistance locus appeared in different 
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Table 4. Resistance loci feature comments. 






Mutation 


n 


% 


Comments form Stanford HIV Drug Resistance Database 


PI Mar - - - 


PI Min 


L10I 


5 


4.7 


L10I/V occur in 5-10% of untreated persons. 




Q58E 


1 


0.9 


Q58E associated with decreased susceptibility to TPV/r and possibly other Pis. 




A71T/V 


16/12 


14.9/11.2 


A71T/V occur in 2-3% of untreated persons. 


NRTI 


M41L 


1 


0.9 


M41L usually occurs with T215Y. confer high-level resistance to AZT and d4T and low level resistance 
to ddl, ABC, and TDF. 




D67H 


2 


1.9 


D67H is a highly unusual mutation at this position. 




D67N 


3 


2.8 


D67N contributes resistance to AZT and d4T. 




T69A/N 


1/2 


0.9/1 .9 


T69N/S/A/I are NRTI-selected mutations. 




K70N 


1 


0.9 


K70Q/N/S/T are rare NRTI-selected mutations. 




K70R 


3 


2.8 


K70R causes intermediate resistance to AZT and low-level resistance to d4T and TDF. 




M184L 


1 


0.9 


M184L is a unusual mutation at this position. 




M184V 


2 


1.9 


M184V/I cause high-level resistance to 3TC and FTC and low-level resistance to ddl and ABC. 




L210W 


1 


0.9 


L210W contributes resistance to each of the NRTIs except 3TC and FTC. It usually occurs with the 
mutations M41L and T215Y. 




T215F/Y 


2/3 


1 .9/2.8 


T215F causes AZT and D4T resistance and reduces susceptibility to ABC, ddl, and TDF particularly in 
combination with M41L and L210W. T215F occurs more commonly (than T215Y) with Type II TAMs. 




K219E 


3 


2.8 


K219Q/E decrease AZT and probably d4T susceptibility when present with K70R or T215Y/F. 


NNRTI 


A98G 


2 


1.9 


A98G reduces NVP and EFV susceptibility by about 5-fold and 3-fold, respectively. 




K101E 


2 


1.9 


K101E causes intermediate resistance to NVP and low-level resistance to EFV, ETR, and RPV. 




K103N 


13 


12.1 


K103N causes high-level resistance to NVP, and EFV. 




VI 081 


2 


1.9 


VI 081 reduces NVP and EFV susceptibility. 




V179D/E 


2/2 


1.9/1.9 


V179D/E occur in —1% of NNRTI-naive individuals and alone reduce NVP and EFV susceptibility. 




Y181C 


9 


8.4 


Y181C causes high-level resistance to NVP, ~2-fold decreased susceptibility to EFV, and ~5-fold 
decreased susceptibility to ETR and RPV. 




G190A 


2 


1.9 


G190A causes high level resistance to NVP and intermediate resistance to EFV. 


doi:1 0.1 371 /journal.pone.0089291 .t004 



clusters. This finding indicates that the prevalence of subtype B 
drug-resistant strains may be the result of both self-replication and 
mutual interaction of viruses in former paid blood donors in 
Henan Province[32-33]. Some PCR sequences with <1% 
difference. This may be the epidemiologic linkage, or strains from 
the same source. Although we tried to obtain detailed information 
from all patients, the information was not always available, 
especially that related to virus spreading. 

The primary resistance rate was 17.7%, which is similar to the 
rate (17.0%) of newly diagnosed individuals reported by Johnson 
et al [20], and higher than the rate (11.6%) reported by Novak 
et al [34], It is also higher than the rate in the United States (10%) 
[35] and the rate reported by Wang X, et al [36] and Liao L, et al 
[37] in China (2.7%-3.8%). However, this rate is lower than that 
in Canada (18.5% in 1999-2000, 27.6% in 200 1) [38-39] . The 
rate of resistance-associated mutations was 40.2%, of which PI had 
no major resistance-associated mutations. The rate of resistance- 
associated mutations of PI minor was 30.8%, and the rate of 
mutations causing potential low resistant was 1.9%. This mutation 
rate is similar to the rate (1.3%) reported by Hamers et al [21], 
and lower than that (1 1.1%) reported by Li et al [40]. The rate of 
NRTI resistance-associated mutation was 10.3%, which is higher 
than the rates reported by Novak (7.8%) and Hamers (2.5%). The 
rate of NNRTI resistance-associated mutation was 18.7%, which 
is higher than the rates reported by Novak (3.0%) and Hamers 
(3.3%) [21,32]. In addition, 6.5% of cases had resistance- 
associated mutations in PI, NRTI and NNRTI, and 8.4% had 



resistance-associated mutations in both NRTI and NNRTI. One 
possible explanation of this phenomenon is that the subtype B 
drug resistance is different from other subtypes, which is consistent 
with the view of Chilton et al [41]. NNRTI resistance was the 
most severe. The rates were in the following order: EFV 14.0% 
(15/107), ETR 10.3% (11/107), NVP 15.4% (16/107), and RPV 
10.3% (1 1/107). NRTI resistance was the next, and the rates were 
in the following order: 3TC 2.8% (3/107), ABC 6.5% (7/107), 
AZT6.5% (7/107), d4T 7.5% (8/107), DDI 6.5% (7/107), FTC 
2.8% (3/107), and TDF5.6% (6/107). PI drug resistance was the 
least severe, with the rates of both NFV and TPV/ r resistance 
being 0.9%. 

The main reasons for the higher rates in the drug resistance 
survey include: 1) People were infected with a single route of 
transmission; 2) The sample size is small; 2) Although the detailed 
interviews were carried out, it still does not rule out the possibility 
that some patients occasionally took antiviral drugs, complain of 
untreated together with treated for a long time. 4) A few patients 
had resistance strains. 5) The high resistance rates may pardy 
result from primary infection by resistant strains, rather than the 
in vivo virus evolution. So transmitted antiretroviral drug resis- 
tance (TDR) in treatment-naive HIV-infected Individuals is very 
important, and needs further in-depth research [42]. 

Characteristics of resistance loci in this study were different from 
the annotations in the Stanford HIV drug resistance database. 
M4 1 L often occurred together with T2 1 5 Y mutation, and caused 
high resistance to AZT and d4T. T2 1 5 Y replaced T2 1 5F, which 
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Table 5. Natural resistance rate of HIV-1 in Henan former paid blood donation population. 









High-level resistance 


Intermediate resistance 


Low-level resistance 


Potential low-level 
resistance 




Drugs 


n 


% 


n % 


n % 


n % 


Protease Inhibitors 


ATV/r 














DRV/r 














FPV/r 














IDV/r 














LPV/r 














NFV 










1 0.9 




SQV/r 














TPV/r 










1 0.9 


Nucleoside RTI 


3TC 


2 


1.9 






1 0.9 




ABC 


1 


0.9 




5 4.7 


1 0.9 




AZT 


4 


3.7 


1 0.9 


2 1.9 






D4T 


4 


3.7 


1 0.9 


2 1.9 


1 0.9 




DDI 


1 


0.9 


2 1.9 


2 1.9 


2 1.9 




FTC 


2 


1.9 






1 0.9 




TDF 






4 3.7 


1 0.9 


1 0.9 


Non-Nucleoside RTI 


EFV 


11 


10.3 


4 3.7 








ETR 






9 8.4 


2 1.9 






NVP 


15 


14.0 






1 0.9 




RPV 






9 8.4 


2 1.9 
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changed intermediate resistance to DDI to low resistance. When 
K103N and VI 081 mutations happened together, they caused 
high resistance to EFV and NVP. Generally resistance loci reduce 
HIV adaptability. But researches revealed some exceptions. For 
example K103N and Y181C mutations almost have no influence 
on the strains, and Y181C mutation even confers the virus better 
replication capability compared to the wild-type strain. Once the 
resistant mutation survives, it will not disappear. Therefore the 
incidence of drug resistance exhibited an increasing trend [43]. 

The mutation rates in the protease and reverse transcriptase 
regions are higher than 60.0%. These results reveal that the 
frequency of HIV mutation is relatively high, and multiple 
quasispecies could exist in the same infected individual. When 
mutations happen in the key regions, a new type or subtype may 
be created. Although the mutation does not directly lead to the 
occurrence of drug resistance, it can reduce the drug-resistant 
genetic barrier, compared with wild type can be selected by 
resistance to faster, when the key drug target nucleotide mutation, 
may lead to a corresponding change in drug sensitivity [44] . 

There are several limitations in this study. Firstly, this is a cross- 
sectional retrospective survey of natural resistance epidemiological 
in former paid blood donors in Henan Province. There was not 
enough information associated with clinical treatment. As a result 
the resistance monitoring results might be biased, and could not 
trace the resistance of individuals. Therefore, it is necessary to 
conduct further follow-up studies for drug treatment effect. 
Secondly, because the sample storage time was too long, only 
about two third of samples could be successfully analyzed. Thirdly, 
because the characteristics (demographic, time, risk factor, 
epidemiological and/or geographic) of the population were 



insufficient, the naturally occurring mutations that clustered 
phylogenetically were not properly described. Fourthly, the high 
prevalence of resistance to NNRTI was not properly discussed 
regarding its main causes and consequences, such as treatment 
coverage of infected population, low genetic barrier of this drug 
class, popular use of NNRTI due to its low cost. Further research 
is needed on the adaptability and replication ability of the 
resistance loci. This may be achieved through follow-up observa- 
tion of the variation trend and the influence of drug-resistance. 
Since the sample size in this study is small, in the future study, we 
will try to increase the sample size. 

In conclusion, through studying HIV drug resistance loci and 
the characteristics of natural resistance in former paid blood 
donors in Henan Province, we understand better the relationship 
between resistance loci and resistance. Understanding the natural 
resistance situation before treatment is important for the 
elucidation of the regular pattern of drug resistance genes, the 
rational use of antiviral drugs, the development of effective anti- 
virus treatments, the selection of salvage therapy, and for 
providing guidance to minimize or delay drug resistance. 
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